subroutine writeSolution(nfile,solution,n)
    dimension solution(n)
    real*8 solution
    do i=1,n
        write(nFile,rec=i) solution (i)
    enddo
end
!*************************************************************************
subroutine readSolution(nfile,solution,n)
    dimension solution(n)
    real*8 solution
    do i=1,n
        read(nFile,rec=i) solution (i)
    enddo
end
!*************************************************************************
!*************************************************************************
!*************************************************************************